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We investigate by the use of the Martin - Siggia - Rose generating functional technique and the 
self - consistent Hartree approximation, the dynamics of the ring homopolymer collapse (swelling) 
following an instantaneous change into a poor (good) solvent conditions. The equation of motion for 
the time dependent monomer - to - monomer correlation function is systematically derived. It is 
argued that for describing of the coarse - graining process (which neglects the capillary instability 
and the coalescence of "pearls" ) the Rouse mode representation is very helpful, so that the resulting 
equations of motion can be simply solved numerically. In the case of the collapse this solution is an- 
alyzed in the framework of the hierarchically crumpled fractal picture, with crumples of successively 
growing scale along the chain. The presented numerical results are in line with the corresponding 
simple scaling argumentation which in particular shows that the characteristic collapse time of a 
segment of length g scales as t* ~ Co<?/ r (where Co is a bare friction coefficient and r is a depth of 
quench). In contrast to the collapse the globule swelling can be seen (in the case that topological 
effects are neglected) as a homogeneous expansion of the globule interior. The swelling of each 
Rouse mode as well as gyration radius R g is discussed. 

PACS numbers: 61.25.Hq,82.35.Lr,36.20.-r 



I. INTRODUCTION 

The equilibrium theory of the coil - globule transition is one of the major issues in polymer physics. Many theories 
have been developed to understand this physical phenomena in more detail. H, |, |, |. The theory for kinetics of 
this transition §, g g |, g, [k| fl2[ |l3| , |l4| , |l5| , [l6| as well as the pertinent experiment |l7j is the subject which has 
recently attracted a broad interest. One of the main motivations for these studies is that the first stage of a protein 
folding process is believed to be a fast collapse to a compact but nonspecific state. This has been shown for example 
by lattice Monte Carlo simulations |ll| . The first stage of folding appears as sequence - independent and therefore the 
process is similar to the collapse of a homopolymer. Generally the collapse of a polymer chain can be addressed to an 
(abrupt) change of the second virial coefficient from the good solvent regime v > to the poor solvent regime v < 0. 
The second virial coefficient depends on the temperature and has a Boyle point v = at the so called 9 temperature. 
The resulting attractive two body interactions requires at least the third virial (always positive) which prevents the 
chain from the collapse to a single point. 

In de Gennes' "expanding sausage model" |6| a flexible chain changes its conformations on a shortest scale through 
the formation of crumples after being quenched to poor solvent conditions. Then crumples are formed on a larger 
scale, resulting in a sausage - like shape, which eventually leads to a compact globule. One can see that this model 
is translational invariant along the chain backbone, provided that the chain has cyclic boundary conditions. This 
simplifies the problem and allows us later on the use of simple Rouse transformations [?, g, of the chain coordinates. 

Concerning de Gennes' model, it is shown by Brownian dynamics simulations |13|, that the "sausages" becomes 
unstable with respect to capillary waves. Consequently a so call ed p earl necklace structure is formed. This brought 
about a number of publications where scaling arguments (To[ [ll| [l2| and the Gaussian self - consistent approach |14f| 
have been put forward. Actually the pearl necklace formation breaks the translational invariance along the chain 
backbone and complicates these issues. 

Alternative recent phenomenological considerations are put forward by joining scaling argument and computer 
simulations j^, |l6) . It was argued |g] that the relaxations times related with the capillary instability as well as with 
coalescence of pearls are short compared to the characteristic time on which the "sausage" change its configuration. 
Since the positions of pearls along the chain are random they can be averaged over. The resulting "sausage" can be 
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seen as an envelop (which has the form of a flexible cylinder) of the pearl necklaces and the overall chain configurations 
are composed by random walks of sausages. This picture successfully explains the so called coarse - graining process 
shortly after the quench to a range of temperature below 8 - point, but above eventual glassy globular relaxation 
regimes ||, This argumentation reconciles in a sense the two scenarios mentioned above, but the microscopic 

picture is still lacking. Therefore instead of considering the processes of formation and coalescence of pearls we 
concentrate here on the coarse - grained dynamics of the "sausage" . This model, as it was already mentioned, is 
homogeneous in the sense that the translational invariance along the chain backbone is effectively assured during the 
process of the collapse. 

In this paper we provide the microscopic theory for the coil to globule transition. In order to get insight and 
some intuitive picture for the later discussion we start from the scaling consideration. Then we study the Langevin 
dynamics of the problem based on the Martin - Siggia - Rose (MSR) generating functional technique together with 
the self - consistent Hartree approximation p2] , p3| , pif . It is quite important for the judgment of the final results 
how the low molecular solvent dynamics is treated. In this paper we restrict our consideration to the random phase 
approximation (RPA) which is well - known in the context of the theory of both polymeric j2fj and low - molecular 
|P6| systems. The RPA fails to account the hydrodynamic interaction because in the hydrodynamic regime collisions 
dominate and keep the solvent in a state of local thermodynamic equilibrium (see e.g. Sec. 6. 5 in ref. pq|), so that we 
leave this subject for the future publications. As a main result the generalized dynamical equation for the collapsed 
(or swelled) chain has been derived. The relaxation laws for the early and late stages are investigated analytically 
whereas the whole numerical solution is also done and thoroughly discussed. 



II. SCALING 
A. Collapse 

Let us first consider the dynamical time scales for the ring polymer collapse under Rouse dynamics conditions. By 
this assumption we neglect certain physical properties, such as capillary instability and long ranged hydrodynamic 
interactions. In this first section, our main goal is to provide a corresponding scaling picture for a Rouse chain which 
we are going to discuss with more refined analytical methods below. Nevertheless the principal times scales are fixed. 
The relaxation of the each length scale is clearly associated with the relaxation of mode as we will show later. 

To model the collapse process we consider the initial state of globule formation to be composed by a Gaussian chain 
of N monomers of size b±. It can be considered as a fractal. In the states at later time the corresponding structure 
is assumed to be the same, generated by a hierarchical random process, as illustrated in Fig.[l]. This self similar 
hierarchical process defines static and dynamic properties by a set of exponents. In each scale and at each stage, 
the chain can be again re-expressed by random walks of N/g coarse grained monomers of size b(g) which contain g 
original monomers |3| . The main effect of the quench to poor solvent conditions can be viewed as if the Gaussian (8 -) 
chain becomes instantly elastic and collapses in order to minimize its contacts with the solvent. At larger times t, the 




FIG. 1: Initial configuration of a chain can be viewed as hierarchical 
fractal in ring geometry. At time t all length scales smaller than a 
certain &2 becomes collapsed. 
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Gaussian fractal structure becomes less complex, since at smaller scales the structure (pi in Fig. [|) becomes collapsed. 
Then all monomers on length scales smaller than, e.g. bi, condense to a scale of length b 2 - On the other hand the 
overall structure on larger scales remains a random walk of the new coarse grained monomers of size &2- After further 
collapse of g ~ b\ monomers, the longest dimension of each collapsed segment remains as 63. The vectorial sum of 
the net force acting on the each segment from outside of the segment is zero if the segment belongs to the fractal 
structure. In the linear chain dynamics, the coarse graining picture based on the fractal structure moves to the late 
stage when the intermediate chain segment can no longer be considered as a part of the fractal due to the continuous 
condensation. The typical late stage configuration for a linear chain is the collinear structure where terminal parts of 
chain experience a net force 

Due to the absence of end effects in the chosen ring geometry, the collapse dynamics is a continuous coarse graining 
process until the polymer reaches its compact globule conformation. We will discuss now an exponent x for the 
characteristic time of the collapse, i.e. ^collapse ~ N x . When the chain is quenched into poor solvent conditions, 
the initial random walk configuration contracts immediately. The energy per each segment of size R\ is E(Ri) ~ 



rgk B T ~ T(Ri/b) 2 k B T where r 



Vf is the depth of quench (yi and V{ are the initial and the final second virial 



coefficients correspondingly). The contraction results in a finite net force, which can be estimated for a given segment 
of length Ri (g) ~ b^fg to be 



/(<?) = 



dE 



knT 



dRi(g) 

The corresponding velocity of each segment is given by 



b 2 



-tRi 



Vs. 
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(2.1) 



(2.2) 



As we assume the Rouse dynamics, the corresponding friction coefficient of each segment scales with the number of 
monomers involved, i.e., £ ~ £n<?, where £n is the bare friction coefficient (determined by the white noise correlation 
and the fluctuation dissipation theorem). Since R\ ~ b 2 g, the characteristic time, t* ~ Ri/u, for the collapse of each 
(coarse grained) segment of length i?i reads 



9 



in 



(2.3) 



Here we define to = Cob 2 /k B T = b 2 / D where D is the diffusion constant of a solvent particle and is related to £n by 
Einstein relation D = k B T/( . Since the chain has self similar structure on all scales, this relation holds until the 
chain merges to a single compact globule. The total time for the collapse is the time required for the largest length 
scale with N monomers. 



N 

f* f 

^collapse ^_ L 0- 



(2.4) 



If we use Zimm dynamics, the friction of the segment of size R\ would no longer be determined by the number 
of monomers, but by the geometric size of the chain |Ts|| , i.e., (z ~ Co-Ri- The characteristic time for collapse is 
t* z ~ i/gto- Note that in Ref.^|, the characteristic time for collapse is t ~ gt$ with Zimm dynamics where the 
infrastructure of a segment is a series of peal necklaces due to the capillary instability. The net force of contraction 
in the geometry of necklace is independent of the scale. The contraction effectively brings monomers which belong to 
the string to the globule. Therefore /necklace ~ k B T /b. We summarize the characteristic time for collapse in different 
regimes in table [I A. 





uniform / ~ ^fg 


necklace / ~ 1 


Rouse C ~ g 




t R ~9 3/2 


Zimm Q ~ yfg 




t n z ~g 



The preceding scenario of the collapse based on the hierarchical fractal picture have been discussed first in ref. || . 
It has been shown there that the theoretical predictions are compatible with the MD - simulation. However we should 
emphasize that at a larger times, when all pearls are merged to a single cluster, the driving energy is no longer 
determined by the fractal regime, E ~ rgksT, but mainly ruled by interfacial effects. This has pointed out by de 
Gennes j(| in his sausage model. The dynamics in this regime is determined approximately by changes of the surface, 
which are described by 



Winter ~ 1 Lr . 



(2.5) 
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In eq. 21 7 ~ fcsT/£ 2 stands for the surface tension, L is the "sausage" length and r is its radius. The sausage is 
formed by thermal blobs of diameter £ ~ b/r. In this regime the collapse can be viewed as a minimization of the 
interfacial energy -Eintcr under the constant volume V ~ Lr 2 . The characteristic time which corresponds to this regime 
is given by |§ 

t** ~ T N 2 t Q (2.6) 

The crossover between the two above-mentioned regimes is detrmined by the match of the two characteristic times, 
i.e., t* fa t** or t 2 N w const. 

This matching is simply understood by taking into account that the globule behaves liquid like under the constraint 
of constant volume. The volume is given by the number of blobs iV/^biob, where gbiob ~ l/ 1 " 2 - Thus we may conclude 
that in the crossover range the number of blobs Nuf^l ~ N/guob corresponds to the equilibrium value. De Gennes' 
dynamics is associated therefore with the re-packing of the incompressible "blob fluid" . 

The presence of these two regimes have been shown by the Brownian dynamics simulation (see Sec. Ill B and Fig. 6 



in ref.[|16|]). It was found that the first regime becomes faster for the larger quenching r (see eq.(2.2)) whereas the 
second one become slow as the quenching grows (see eq.( [2.6| )). 



B. Swelling 

When the solvent quality is changed suddenly from poor to good or 9 - conditions the globule starts to swell. We 
consider here for brevity a quench from a poor solvent to 6 - solvent conditions and assume that the globule interior 
always contains enough solvent molecules i.e., the solvent is not very poor (r < 0) . Then one can ignore a specific role 
of the solvent transport and assure the homogegeneous globule expansion. In this subsection we restrict ourselves to a 
scaling picture to set up the main time scales by employing an "expanding blob" picture. Let us therefore assume that 
during the process of expansion the segments on the length of blob size £ maintain local equilibrium states, whereas 
on a larger scales the system is no longer in equilibrium. This amounts to the assumption that the initial state of the 
globule can be described by free energy which is proportional to the number of blobs times the thermal energy [0. 
We can then write the free energy for the overall quasi - equilibrium regime in the following form 

^giob oc k B T (-) , (2.7) 
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where the number of monomers in the blob, g, is a time dependent function and still has a Gaussian statistics, i.e. 

• < 2 - 8 > 

We assume as well that the quasi - equilibrium scaling for the overall globule size, R, is valid and gives 

/iV\ 1/3 

Rozd-jj . (2.9) 



The substitution of eqs.([2.g|) and (UM in eq. (2.7) yields 



b 6 N 3 

Fgiob <x fc B T— . (2.10) 

To take into account the correct long time limiting behavior, when the chain size becomes proportional to b 2 N , the 
elastic energy term should be included. This term counterbalances the over-swelling and stabilizes the system, so that 
the whole quasi - equilibrium free energy reads 

/b 6 N 3 R 2 \ , , 

F «** T {-w + m) • (2 - n) 

After that the corresponding equation of motion for R(t), i.e. 

C^ = -^ (2.12) 
^ dt SR y ' 
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takes the following form 



N dR(t) _ b 6 N 3 



D dt 



R 



R 7 b 2 N 



(2.13) 



In eq.(2.12) £ ~ (qN is the Rouse friction coefficient and in eq.(2.13) D = ksT/(o denotes the monomer diffusion 
constant. It should be noted, that we do not use the Rouse diffusion constant throughout the paper, since we wish to 
keep track of the chai n len gth dependences explicitly. 



The solution of eq. fl2.13|) is 



R 2 (t) 



p8 -8Dt/b 2 N 2 

tt e 



+ 6 8 JV 4 ( 1 - e SDt/b 2 N- 



1/4 



(2.14) 



In the asymptotic limit of t — > oo, we recover the Gaussian chain size, R 2 ^ — b 2 N . The characteristic time for swelling 



is determined by the condition: R 



1/2 



b 2 N . The characteristic time t* 
f* ~ N 2 

''enroll 1 " 



ell scales with the chain length as 

(2.15) 



We will see in Sec.V, devoted to the numerical investigation of the full equation of motion, that because of coupling 
of different Rouse modes the driving force of the swelling is mainly determined by the large length scale (see e.g. 

In this consideration the chain is treated as "phantom" , notwithstanding the effective interaction is taken into 
account through the virial coefficients segments still can cross each other. This use of pseudo potentials does not allow 
to include topological effects discussed recently. If, however, the effects of topological constraints are dominant, the 
chain can be swollen only by reptation through the channels made of neighboring segments. The swelling process of 
an entangled globule needs more consideration J27], ^ . 



C. Effective Hamiltonians and local time scales 



The scaling shows what can be expected at large length scales only. It is at first sight difficult to estimate the 
expected time scales on the level of Rouse type modes. Nevertheless some remarks to the scale dependence on the 
relaxation modes can be drawn. In the more detailed theory, which is presented below one can recover the dynamical 
scales on the very late relaxation regimes. The chain size in globular states scales as R ~ &(iV/r) 1 / d for scales larger 
than the thermal blob size £t = b/r. A general staring point is the Edwards Hamiltonian H(TL(s)), where R(s) defines 
the chain position vectors and s the contour variable. Any dynamic theory can be formulated through the Rouse modes 
which are usually defined by a Fourier transform of the position vectors as X(p) = (1/N) R(s) exp(isp), where p 
stands for the Rouse modes with p = 2nn/N and n = 0, 1, 2, . . . , JV — 1 |lj|. In general the effective Hamiltonian for 
a collapsed chain can be written in terms of Rouse modes as a simple quadratic form 



H cS = k B Tf( T )Yy a \X( P r , (2.16) 



where a is determined by the static correlation and the function /(r) is 1 for the good and 9 - solvents but /(r) ~ T 2 l d 
for poor solvent conditions. In the case of extended chains it is easy to show that 2a = 1 + 2v ( where v = 3/(d+ 2)) 
whereas in the case of compact chains 2a = 1+2/d. The Rouse dynamics for these effective chain variables X(p) is 
determined by an Rouse equation of the form 

Co ^M + !^ /(r) ^X (p) = (2.17) 

This effective Rouse equation set the time scales for the latest stages of the dynamical evolution. Namely, for the 
swelling the relaxation spectrum has the form 



whereas for the collapse it is 



T "*(P) * T 2/ d % Vd ■ ( 2 - 19 ) 



This corresponds naturally to the characteristic Rouse relaxation times as we will discuss below in more refined 
theories it in Sec. IVB. 
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III. EQUATION OF MOTIONS FOR CORRELATORS 

A. Model 

In this section we provide a more general formulation of the Langevin dynamics for a polymer chain immersed in 
the solvent. The chain conformation is characterized by the d - dimensional vector - function R(s, t) of the time t 
and s (1 < s < N), the segment position along the chain contour. The corresponding intra chain Hamiltonian has 
the following form 



1 N 

H=^J2^sR( Sl t)} 2 +H int {R( S ,t)} , 



(3.1) 



where e = dkgT/b 2 is the elastic modulus with the Kuhn segment length b, N is the general number of segments and 
the finite difference V s Rj(s,<) = Rj(s + l,t) — R,j(s,i) and the intra-chain interaction Hamiltonian has the form 



N N 



H int {R(s,t)} = i£)5S(R( a ,f)-R(a',t)) 

s=l s'=l 
N N N 



(3.2) 



s=l s' = l s'' 



In eq.(3.2) v(r) and w(r 1; r 2 ) are the second and the third virial coefficients correspondingly. Let us treat the low - 
molecular solvent molecules as a separate component and specify their local positions by the vector - function r( p ) (t) , 
where p = 1,2, ... M numerates the number of the solvent molecules. We denote by V ps (r) and 14s(r) the polymer - 
solvent and solvent - solvent interaction potential correspondingly. After that the whole polymer - solvent dynamics 
is described by the following Langevin equations: 



d 

Co-^Rj(s,t) - eA s R j (s,t) + 



SRj(s,t) 



flint {R(S,*)} 



5> ps (R(M)-r«(t)) =/i(*.*) 



(3.3) 



and 



C— r {p) (t) 



N 



+ 



J2 V ss (rW(i)-rW(^ 

71=1 

N 

J2V PS (rW(t)-R(s,tj)=fM(s,t) 



(3.4) 



where j numerates Cartesian components, Co and ^ are bare friction coefficients of polymer segments and solvent 
molecules (which should be of the same order and are put equal in the following) and the second order finite difference 
A a Rj(s,t) = Rj(s + l,i) + Rj(s- l,t) - 2Rj[s,t\. 

In order to reformulate the Langevin problem (p^) - ( ]3 .4| ) in a more convenient form we go to the MSR - functional 
integral representation fl24|. The generating functional (GF) of the problem has the form 



Z{---} = j DRj(s,t)DRj(s,t) 
x exp|s R(s, t),R(s, t) 



R( s ,t),R( s ,<)]} , 



(3.5) 
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where the influence functional (which describes the influence of solvent molecules on the chain) reads 

In f JjDrW(t)i»W(t)exp J Advent [r {p \M 
J P =i { 

UN - 

6 



p— 1 s— 1 

A-/ iV 



£E / (>)(*) - R(M)) 

— i — i ^ OTj \t) 



p—1 s—1 ' 



the intra-chain action is given by 



R(s,t),R(s,t) = E / (ftsi-Rj(M) 



Co^i(s,i)-eA s ^(s,f) 



5Rj(s,t) 



and the solvent action has the form 



M 



.4. 



solvent 



rW(t),rW(t)] =Ey*| if ? ) (*) 
+ fc fl TCoH p) (t) 



i*r?°(') 



-r 



(3.6) 



(3.7) 



(3.8) 



The representation (3.5) - ( |3.8| ) is a suitable starting point for further approximations. In particular these expressions 
are very convenient for integration over collective solvent variables, which will leads to the effective "actions" for the 
polymer. We remind that this procedure is only possible if the solvent particles dynamics is much faster compared to 
the polymer dynamics. While this is in general ensured for long chains, nevertheless this point has to be treated with 
a care for the polymer collapse problem. 



B. Self - consistent Hartree approximation 



As mentioned already it is an important to point out how the solvent dynamics is treated. Indeed we may follow 
at least two different ways. The solvent can be considered within the hydrodynamic approximation in the same way 
as it was done in ref. |29j| , where the solvent dynamics was described by an incompressible Navier - Stokes liquid. 
This approach leads to a time - dependent hydrodynamics interaction. The other way is to treat the solvent as a 
dynamical background in a random phase approximation (RPA). Here we restrict ourselves to the dynamical RPA 
pi] , |30| which is well known in low molecular liquid dynamics . It is also known |2(| that the RPA is a mean field 
type description in which the free - particle behavior is modified by an effective interaction. The main drawback of 
this approximation is that it does not give the proper hydrodynamic behavior. The reason for this lies in the fact 
that RPA neglects collisions which dominate in the hydrodynamic regime. We will come to this point in a future 
publication. 

In order to accomplish the calculation for the present purpose we make use of the the transformation to the collective 
solvent density |2^| and integrate them out. The level of this procedure will determine the level of approximation 
also. We relegate all technical details of this calculation in the Appendix A. 

The Hartree approximation will then take into account all mean field diagrams. Naturally the mean field description 
appears poor in good solvent conditions, but becomes better in the globule state, since fluctuations are less important 
under the increasing globule density. 

The GF which is determined by eqs.([A5|) and ( |A6| ) is still highly nonlinear with respect to R(s,t) and R(s, t). In 
order to handle the difficulties which are associated to these we use Hartree - type approximation. In this approxima- 
tion the real MSR - action is replaced by the Gaussian one in such a way that all terms which include more than two 
fields R(s, t) and R(s, t) ) are written in all possible ways as products of pairs of R(s, t) or/and R(s, t) coupled to the 
self - consistent averages of the remaining fields. In ref. it was shown that if the number of field components is 
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large the Hartree approximation and the next to the saddle point approximation merge and both become exact. The 
resulting Hartree action is a Gaussian functional with coefficients which could be represented in terms of correlation 
and response functions. All calculations are straightforward and details can be found in the Appendix B of ref. p|. 
The only difference is that here the second and third virial terms (two last terms in eq.(A6)) explicitly enter into the 
equation. After the collection all terms the final GF reads 



£{•■•} = J DRDR cxp {4°L [ R > H-] 



N N roc 



EE 

s=l s'=l 
N N 

EE 

s=l s' = l 
2 p°° 



dt 

— oo J — oo 



dt 



dt 



dt 



dt' iRj (s, t)Rj (s\ t')X(s, s'; t, t') 

dt' iRj(s, t)Rj(s', t)X(s, s'; t, t') 
'iRjis^iRjis'^x&s'it^')} , 



where 



X(s,s';t,t') = ~G(s,s';t,t') J ^^k 4 \V ps (k)\ 2 F(k; S ,s';t,t')Soo(^t,t') 



d d k 
(2tt)° 



k 2 [|^ ps (k)| 2 ^ i(k; t, - V pa (H)6(t - t')] F(k; s, s'; t, t') 



and 



N 

E 

s" = l 



X(s,s';t,t') 



d d kd d q 
(27r) 2d 



w(k, q)F(q; s' , s"; t, t)F(k; s, s'; t, t)6(t - t') 



d d k 
(2^ 



^fc 2 |T/(k)| 2 F(k; S , S ';t,t')5oo(k;i,i') 



In eqs.( p.9[ ) - (3.11) the response function 

G{s,s';t,t') = (iB,{s',t')K{s,t)' 
and the chain density correlator 

F{k;s,s';t,t') = exp |-^-Q( s , s '; t, 



with 



Q(s,s';t,t') = (R{s,t)R(s,t)) - (R(s, t)R(s', t')) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



In eqs.( [3.9| ) - (3.11) 5po( k; t, t') and Soi(k;t,t') are the solvent RPA - density correlation and response functions 
correspondingly (see eqs.(A7) and (A8) in the Appendix A). They embody information on the s olve nt dynamics. 

The pointed brackets denote the self - consistent averaging with the Hartree - type GF, eq.(3.£}. Below we will 
also concern transient time regimes, so that keeping both time arguments for correlator F in eq.(3.1C) equal to each 
other does not necessarily mean that this is a static correlator F st . On the other hand we assume that the fluctuation 
- dissipation theorem (FDT) hold for both chain and solvent correlators, then 



d 

G{s,s'\t-t l ) = {k B T)- 1 — Q{s,s'-t-t') at t > t' 



(3.15) 



Soi(k;t-t') = (k B T)- 1 -^S o{k;t-t') at t > t' 



(3.16) 



Note that in eq.( 3.16| ) the units of the correlation function 5oo an d the response function Sqi are different (see the 
corresponding Appendix A for the notation). 
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t', we obtain 



Z{-..} 



Now we can use eqs.(3.15) and (3.16) in eqs.(|3.9|) - (3.14). After integration by parts with respect to time argument 



DRDKexp 



N 
s.s' — l 



(It 

oo J — oo 



d 

dt' iRj (s, t) [C S(t - t') + 9(t - t')T(s, s'; t, t')} —Rj (s\ t') 



N 



oo f-t 

^2/ dt dt' iRj(s,t)n(s,s';t) Rj(s',t) 

s ,s' = l"'- 00 "'- 00 

N pOO poo 

k B T / dt / dt'ii? i (s,t)[Co<5(i-i') + 6 '( i - i/ ) r ( s ; s/ ; t > i ')]^j( s/ ^')} 



where the memory function 



T(s,s';t) 



d d k 



k B T J (2tt) 



fc 2 |y ps (k)| 2 F(k; S , S ';^oo(k,t) 



(3.17) 



(3.18) 



and the effective elastic susceptibility 

n( S) s';t) = e6 sa ,A s -J ^ fc 2 V(k) 



N 



F(k; s, s';t, t) - <W ^ F(k; s, s"; t, t) 



s" = l 



1 N 

E 



{2n) 2d 



k 2 w(k, q) 



F(k; s, s'; t, t)F(qj s", a'; t, i) - <W ^ F(k; s, s'"; t, t)F(q; a'", s"; t, i) 



In eq.( |3.19 ) the Fourier components of the effective segment - segment self - interaction function is given by 

|l/ ps (k)| 2 $ st (k)/fc s T 



V(k) = v(k) 



i + y ss (k)$ 8t (k)/fc B T ' 



(3.19) 



(3.20) 



where $ s t(k) is the static density correlator for the free solvent s ystem . In eq.( 3.20 ) the second term results from 
a coupling with solvent degree of freedom. The memory function (3.18) describes the renormalization of the Stokes 
friction coefficient £ which originates from a coupling between polymeric and solvent fluctuations. The effective 
elastic susceptibility, eq.(3.19), gives an accou nt of t he non - dissipative contributions, which arises not only from the 
local spring - interaction (the first term in eq.(3.19)) but also from effective two and three point interactions. 



C. Equation of motion 

Now we are in a position to derive equation of motion for the time - displaced correlator 

C(s,s';t,t') = (R(s,t)TL(s',t')) 

as well as for the equal - time correlator 

P(s,s';t) = (R(s,t)R(s',t)) . 



(3.21) 



(3.22) 



The standard way to derive equations of motion by starting fr om the Hartree action (see eq.(3.9) is discussed in 
Appendix B ref.j22|. Using this way for the Hartree action in eq.(3.9) and taking into account FDT yields 



d 



N 



Co—C(s,s';t,t') - J2 n{s,m;t)C(m,s';t,t') 



m—1 
N ,t 



r d 
V / T(s,m-,t,T)—C(m,s';T,t')dT = -2k B T( G(s',s;t',t) 



(3.23) 
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In the case of t' < t the r.h.s. of eq.( 3.23j ) is zero and the time - displaced correlator satisfies the equation 



d N N r* 

(o—C{s )S ';t,l/)-J2 ft( S) m;t)C(m,s';t,l/) + J2 / r(s,m;t, 

m=l ra=\^ 1 ' 



T) — C(m,s';T, t')dr = 

OT 



In order to derive the equation for the equal - time correlator ( 3.2 2| ) we recall that 



g t P(s,s>;t) 



-C(s,s';t,t') 



g^C(s,s';t,t') 



and the initial condition 



CoG(s, s'; t + 0+ t') = -dd(s - s') 



(3.24) 



(3.25) 



(3.26) 



Let us make the permutation of time moments, t^_t', in cq. (3.23). Combining this equation with the original one 
and using eqs. (3.25) - ( 3. 26] ) in the limit i = t' + eate^0, one can derive the result: 



N 



(3.27) 



m — 1 



It is of interest that in eq.(3.27) the memory term is dropped out. 

As it was discussed in the Introduction , here we restrict ourselves to the case where the translational invariance 
along the chain backbone holds during the collapse (swelling). Generally speaking, the presence of the "pearl necklace" 
structure breaks down this invariance, so that the correlator P(s,s';i) depends not only on the "chemical distance" 
\s — s'\ but also on the position along the chain backbone. Even the interface might violate this invariance because 
the chain segments on the surface and in the bulk experience quite different environment. Nevertheless, in the case 
when the pearl formation is fast compared to the dynamics of the envelop and there positions along the chain are 
random these differences can be averaged out (by preparing an appropriate ensemble of pearls realizations) and the 
effective invariance still holds. In this situation the Rouse components are the "good variables" and it is worthwhile 
to make the Rouse transformation in the standard way K9|: 



1 N 



(3.28) 



and 



271 



p=0 



(3.29) 



where the Rouse mode p — 2ttti/ N at n — 0, 1, . . . , iV — 1 and we have used for simplicity the cyclic boundary 
conditions. After that the eq.(3.27) reads 



(2D)- 1 ^- t P(p;t) + n(p;t)P(p;t)^dN- 1 , 



where D — kgT/CQ is the bare diffusion coefficient and 



~/ 2d, N 

tltot) = ^(l-cosp) + ^ 



N 2 



d d kd d q 



fc 2 w(k,q) F{k;p;t,t)F(q;p = 0;t,t) 



2k B T J {2w) 2 d 
- F(k;p = 0;t,t)F(q;p = 0;t,t) . 

On the other hand, the Rouse transformation of the chain density correlator has the form 



1 N f 
F(k; p; t, t) = — 2J cos(pn) exp <^ 

n=l 



-Q(n;t,t) 



(3.30) 



(3.31) 



(3.32) 
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For the short range segment - segment int eract ion one can neglect the wave vector dependence in V and w(k, q) by 
putting V fts v and w(k, q) w w. Using eq.( 3.32| ) in eq.(3.31) and performing the integration over k and q yields 



2d, 



. d2+ 2 1 — cos(pn) 

fi(p;*) = TJ 1 -cosp -f > — 

& 2 2& B T(47r)2 ^ [Q(„, t )]-i- 



4fc B T(47r)' 



N N-n 

EE 



1 — cos(pn) 



[Q(n,t)]^[Q(M)P 



(3.33) 



We stress that the equation of motion (3.3C), ( 3.33| ) for P(p;t) is highly nonlinear because the correlator Q(n,t) 
depends on all other P{k\ t) (which provides also a Rouse mode coupling) as follows 



Q(n,t) = P(n,n;t) - P(n,0;t) 

[1 — cos(Kn)] P(k; t) 



2iT 

E 

k=2tt/N 



In equilibrium all equal - time correlat ors in eqs.( 3.33| ) - ( [3.34| ) can be seen as static ones: P(p;t) 
Q(n,t) — ► Q st (p). Using this limit in eq.(3.30) leads to the following equilibrium equation 



(3.34) 
C st (p) and 



[NC Bt (p)Y 



b 2 



— cosp) 
d d+i 



^ 1 



cos(pn) 



2fc B T(47r) 

JV N-n 



/ j d + 2 

n=l [Qst(n)] 2 



d E E 



1 — cos(pn) 



AkB T^)d^^ [Q st (n)]^ [Q rt (n)]* 



(3.35) 



This equation is identical (with the accuracy of prefactors) with the variational (Euler) equation which we have 
recently discussed in ref. J33, 34 1. This provides the means for answering one of the important questions of whether the 
Edwards Hamiltonian can be equally used for dynamical calculations. Within the Hartree approximation the answer 
is positive provided that the second and third virial coefficients are considered as free parameters. In fact, the Hartree 
approximation is a dynamical counterpart of the variational approach 1 33 , p4[ . 



IV. EARLY AND LATEST STAGES OF THE COLLAPSE (SWELLING) 

Now we are in a position to consider some limiting cases which allow an analytical investigation. These arc obviously 
very early and latest stages of the transformation. We shall defer the full discussion which based on the equation of 
motion numerical solution until the next section. 

Let us consider the abrupt solvent quality changing between the initial v\ > and the final Vf < 0, which is the 
case of the collapse experiment. For the swelling case these are the following: Uj < and Vf > 0. 



A. Early stages 



For the very early stage the eq.( 3.30 ) can be linearized around the starting state, which leads to the following form 

(2D)- 1 ^- t P(p;t) + n st (p)P(p;t)=dN- 1 , (4.1) 

where 

2d 



d2+ 2 1 — cos(pn) 
— (1 - COSJ3 - Vf d > — 



fist (?) 



d d+l 



N N-n 



EE 



1 — cos(pn) 



AJsbTW^^ [Q rt („)]^ [Q st (n)]i 



(4.2) 



In the eq.(4.2) the static correlator Q s t(n) should be seen as a coil correlator, i.e. Q s t(n) = Qs° ll ( n ) m the case of 
collapse and as a globule correlator, Q s t(n) = Qst° b { n )i m the case of swelling. 
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Combining eq.(4.2) with eq.( |3.35 ), which is valid in the equilibrium at v = Vi, yields 



n st ( P ) = d[NC st ( P )}- 1 + ( V - i -v ( )-^-± l-oo S (pn) 
The very early stage can be characterized by the initial relaxation rate 



X(p)=^P(p;t) 



= [dN- 1 - n st (p)C st (p)] (2D) , 



t=o 



which with the use of eq. (^Ij) can be represented in the following form 

(Vi - Vf) 



Co 



Cstip) Hp) 



where 



rf2+ 2 l-cos(pn) 

^0) = 7TTT 

47r 2 n =i 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



The relaxation law for the decrement of the gyration radius ARg(t) = J2 p =2it/n [P(P'-> ~ P(p\ 0)] a t the early stage 
takes the form 



AR 2 (t) 



(vj - Vi) 
Co 



[ E c ^(p)Hp) 

p=2ir/N 



(4.7) 



Taking into account the asymptotic behavior, [C s t(p)] 1 oc p 1+2u (for small p), Q s t(n) cx n 2v (for large n), where v is 
the Flory exponent, one can obtain for function T(p) the following scaling 



(4.8) 



For the collapse case U; > Vf. For quenching from the good solvent, v = 3/(d + 2), and T c (j>) cx p 2 , whereas for the 
quenching from the - solvent v = 1/2 and T c ij)) cx p d / 2 . 

In the case of swelling wj < «f. If one heat up the system starting from the globule state then v = 1/d and 
F s {p) oc p 2 ^ d , whereas for the 8 - solvent initial state one has ^F s (p) oc p d / 2 . 

We will leave the discussion of the relaxation rate X(j>) till the next section but one can immediately see that X(p) 
has the following scaling forms. For collapse 



X c (p) oc 



Vi-Vf J l/ p (4-d)/(d+2) 

\/ P 2 - d / 2 



good — solvent 
8 — solvent 



(4.9) 



where the first and the second lines refer to the quenching from the good and 8 - solvents correspondingly. 
For the swelling process we find respectively for different initial starting points 



X s (p) 



Vf - V; 



1/P 



, poor — solvent 
Co 1 l/p 2 ^/ 2 , 6- solvent 



(4.10) 



Here the first and the second lines are assigned to the heating up from the poor and 8 - solvents correspondingly. 



B. Latest stages 



Now the system is close to the g lobu le state (if the collapse is under discussion) and we can linearize the equation of 
motion around it. In this case eq.(4.1) is still valid but f2 s t(p) should be calculated in the globule state, which comes 



out of eq.(4.2) after substitution Q s t( n ) ~ 
eq.C||) at v = v u C st (p) = C^ l) (p) and 



Q^ l '(n). I n its turn correlator for the globule state Q^ l> (n) satisfies the 



tin) 



Qifin). 



After simple calculations one can obtain 



nif( P ) = d NC^ip) 



(4.11) 
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where 



NCW (p)j oc p 
= — i?o(oo) we have 



Q n i a ^ es ^ stage only the relaxation mode p = 2n/N contributes. Then for 



Ar 2 (t) ocexp<^ -2D fl 



2tt , 

TV 



(4.12) 



i.e. the characteristic relaxation time is given by 



r rcl oc ±N 1+ V* 



(4.13) 



which agrees with ref. M. Note that we derived this behavior as well previously when we had used the effective 
Hamiltonian for the collapsed chain. Obviously these crude estimated describe already the main issues of the dynamic 
behavior. 

The same remark holds for the case of the large time limit for the swelling in a good solvent (vf > 0). The formal 
dynamic theory describes these mode dependence as 



1 d + 8 

Trel OC — 



(4.14) 



which is also consistent with ref.||. Again we derived this already by the effective chain Hamiltonian in Sec. IIC 
Physically this appears naturally since the choice of th e exponen t in e q.( 2.16| ) uses implicitly the FDT to determine 
the statics correctly. It is interesting that eqs.(4.13) and eqs. Q4.14 ) can be seen as a characteris tic R ouse time, 
TRouse ~ N 1+ 2v (see e.g.[p5|), where the Flory exponent v = 1/d or v = 3/(d + 2) in the case of eqs.( |l.l3] ) (collapse) 
and eqs.(L14) (swelling in a good solvent) correspondingly. 



V. NUMERICAL STUDIES 



A. Collapse 

To analyze these results in more detail we performed numerical studies of the corresponding equations. We studied 
so far only the asymptotic behavior, w herea s the intermediate time regimes are not accessible. To do so we explicitly 



computed the numeric solution of Eq. 3.3C for chain length N = 2 6 — 2 10 . For the whole procedure, the three body 
interaction term is fixed as V3 = b e = 1. The time scale is measured in units of t$ = ^b 2 /ksT. The monomer density in 
a globule is determined by the balance between the two-body attraction and three body repulsion, p = \v2 |/i>3 = —r/b . 
The excluded volume provides the condition for the maximum density /o max |i'2| < 1- 

As mentioned earlier, the chain is phantom, i.e., we use pseudo-potentials only which might cause a faster collapse 
since the chain segments are allowed to pass through each other. However, this effect is maybe a minor correction at 
least in the early stages at the beginning of the collapse. At later stages artifact caused by the phantom assumption 
(such as an overshooting on the relaxation curves) are expected to be more pronounced. In order to assure the validity 
of the virial expansion and to minimize these artifacts, the solvent quality should remain in the limit where i>| /t>3 <C 1. 
Let us now discuss the results in more detail. 



1. Solvent quality dependence 

First we are going to present the dependence on the solvent quality as a first check of the theory. The mean squared 
radius of gyration R 2 (t) for various second virial coefficient v-i is plotted in Fig. || for TV = 512. Due to the finite size 

effect the solvent quality which gives the coil to globule transition is expected to be shifted by ~ \j\f~N . We observe 
the transition occurs around r k, —0.25. In equilibrium, the size of the globule in the poor solvent regime varies as 
R ~ b(A^/|-r| ) 1 / 3 , according to the solvent quality. This is in good agreement with the scaling predictions. 



2. Chain length dependence 



The next issue is the chain length dependence. In Fig. ^, the mean square radius of gyration R 2 {t) during the 
collapse is shown for various chain size N. The cha ract eristic total collapse time increases linearly as chain length N, 
as predicted from the scaling estimate given by eq.(2.4). 
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FIG. 2: The mean square radius of gyration Rg(t) as a function of time after quench to poor solvent conditions from 8 solvent 
condition where the chain configuration is Gaussian. The equilibrium R 2 s {t) is shown in inset for various solvent quality. The 
coil to globule transition is observed at solvent quality |r| = 0.25. 
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FIG. 3: The mean square Rg(t) of various chain length N as a function of time after quench to poor solvent conditions 
(r = —0.4) N = 64, 128, 256, 512, 1024 from bottom to top. The inset shows the total characteristic time for collapse for each 
chain vs. their length N. 



3. Relaxation times for individual modes 



During the coil - globule transition the collapse occurs in hierarchical patterns, such that the collapsed segment 
on a smaller length scale can be considered as the unit length of the larger scale of the chain. In the process of the 
collapse transition the active modes disappear one after another. We expect such a hierarchical collapse to be visible 
in the corresponding Fourier mode correlator P(p, t) where p — 2im/N and n — 1,2,..., iV — 1. The characteristic 
time for mode p is related to the relaxation time of the subchain of length g = N/n After the elapse of time t n all 
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length scales less than g are collapsed (see the Sec. II) . There exist N / g collapsed sub-chains at a given time and the 
hierarchically larger length structure is a random walk of these sub-chains. At times less than t n , the contractions on 
length scale smaller than g contribute essentially to the decrease of the correlatator P(p ~ l/g,t). For times t > t n , 
further contractions from the larger scales are reflected. Indeed, the relaxation behavior of P(j>, t) on smaller length 
scales (or larger n > 2) shows two clearly distinctive regimes, see Fig. [|. However, the longest mode (n = 1) shows 
only a single slope until the chain finds itself in the globular state. We may define the fast decreasing regime as the 
first characteristic time scale t* for each mode p. The first characteristic time t* is therefore related to the "internal" 
relaxation time of subchain of length g = N/n. The subchain relaxation time scales with its length g as t* ~ g ~ 1/p 
which is consistent with the scaling results (see Sec.IIA) . The size of each length scale continue to decrease after time 
t* and this is due to the larger scale contraction. 



It is now instructive to give a more general estimate for the characteristic time t*. According to eq.(4.5) the driving 



force for the collapse transition scales as f p ~ (u; — V{)C st F(p) ~ rp' 



vd-2 



t/9 



ud-2 



In d = 3 and at v 



fp ~ T ^/9~^ which agrees with eq. (|2.lD . In the same manner as in Sec.IIA, the characteristic time scales as t* 



where R describes the chain properties as R 



and u ~ fp/Cog, so that the resulting scaling reads, 



1/2, 
R/u 



t* o,^ ff "(i+<*)-i 

T 



(5.1) 



For three dimensions, d = 3 and v — 1/2 we recover the relation, t* ~ Co9/t, which is supported by the scaling 
analysis in the Sec.IIA. 

After the time t*, each Fourier mode decreases slowly (for n > 2). The total collapse time for each mode is identical 
with that of the first characteristic time for the longest mode, i* =27r /jv ■ The slope decrease after time t* is controlled 
(because of modes coupling) by the contraction on longer scales (or smaller p modes) . The fast relaxation time for 
the second largest mode for the chain length N coincides with the total collapse time for a chain of the corresponding 
half length (see again the Sec. II). Naturally, the longest relaxation mode has only one characteristic regime. 
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FIG. 4: The relaxation of each mode P(p = 2irn/N,t) during the collapse transition. 



B. Swelling 

We turn now to the the case of swelling and present the numerical results of the globule expanding after a quench 
from poor to a good (or 9 - solvent) condition. The initial configuration is now assumed to be a compact globule which 
is prepared by a poor solvent with a certain negative second virial coefficient. Contrary to the collapse dynamics, 
which based on the hierarchically crumpled fractal picture, the swelling can be conceived as an homogeneous extension 
of the different Rouse modes. 
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The time dependence of P(p, t) is shown in Fig. [| for different values of the mode index p. Small length scales 
(or large p) relax fast, in the sense that they are saturated faster, while larger scales grow slower. There is no 
clear dynamic exponent observed for the subchain relaxation time in the numerical solution. In the beginning of the 
swelling, overshooting is observed for large p values. This is an obvious artifact of a phantom chain. The monomers 
move out from the dense globule explosively as soon as the solvent condition is switched to a theta solvent. The 
phantom chain can pass through surrounding dense phase without being hindered by the topological constraint. The 
overshooting shows up when the size of the subchain of g monomers r(g) is comparable to the boundary size of the 
total chain R(t). It appears also that the virial expansion is not simply applicable for at least the early stage of 
swelling. 

10 5 
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ci 
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10 3 
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10000 20000 30000 40000 50000 

t/t 

FIG. 5: The mode correlation function P(p,t) for each mode pN/2n=l, 2,4, 8, 16 from top to bottom. The solvent quality is 
switched from vi = -0.3& 3 to v t = -0.16 3 . The chain length is N=1024. 

The growth of the overall size B? g (t) ~ t z is also computed. The dynamic exponent is approximately z ~ 1/2. The 
total relaxation time t* wel , grows as iV 2 as the chain length increases. 

In order to demonstrate how different modes contribute to the swelling, it is useful to calculate the normalized 
relaxation rate X(p, t) = 1 — P(p, t)Sl(p, t)/d as a function of wave vector p. We have discussed already (see Sec. IV) 
the asymptotic forms, i.e. at p <C 1, of this function at the initial time moment. In Fig. |?], we compute the relaxation 
rate X(j>, t) for the swelling condition. When the system is suddenly quenched to the 9 solvent, the compact globule 
conformation is unstable and far from the equilibrium. The relaxation rate X(p,t) reflects the discrepancy between 
the desired conformation and the conformation at given time t. If the system tends to the equilibrium, the relaxation 
rate obviously vanishes. The contribution to the total swelling from different modes p changes as swelling proceeds. 
The length scales smaller then the blob size do not contribute in the swelling. We may define p max = 27r/<? as the 
Rouse index of the shortest active mode at t = 0, where g is again the number of chain units in the Gaussian blob. 
(The larger p > p max do not contribute). As the chain swells, the size of the boundary is growing, but on the other 
hand, the position of p ma , x (t) decreases accordingly. In Fig. [7] one can see the shifting of p ma , x (t) to the smaller value 
of p in the course of time. These findings are consistent with the "expanding blob" picture which we discussed in the 
Sec.II B. 
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FIG. 6: R g (t) for chain length N = 256, 512, 1024. 
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FIG. 7: The relaxation rate 1 — Q(p)P(p)/d for half cycle of p with chain length N=1024. The contribution to the swelling 
comes from each mode p < p c at given time. The positive value regimes are active mode. The inset shows the log-log plot of 
the relaxation rate. 



VI. CONCLUSION 



Using the MSR - generating functional technique and the self - consistent Hartree approximation we have derived 
the equation of motion for the time dependent monomer - to - monomer correlation function. The numerical solutions 
of this equation for the Rouse modes and the gyration radius R g in the cases of collapse and swelling regimes have been 
discussed. It has been shown that the simple scaling arguments for the collapse based on the hierarchically crumpled 
fractal picture match pretty well with our numerical findings. The size of the crumples (or the size of the collapsed 
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segment of the chain) grows successively, so that the characteristic time of the collapse changes linearly with the 
length of the collapsed segment. Swelling, on the contrary, goes homogeneously, with the Rouse modes relaxing with 
different relaxation rates. The spectrum of the relaxation rate spans the interval between the minimal p m i n = 2ir/N 
and maximal p max = 2n/g Rouse mode indices, where the "expanding blob" length, g, is a growing function of time. 

In this paper we have restricted ourselves to the ring polymer. The presence of free ends in the open polymer 
chain brings some special features mainly on the later stage of the collapse. This question have been discussed in the 
number of papers ^ pq |. In particular in ref. || it has been shown that the late stage configuration consists of two 
globules with connecting bridge between them (dumb - bell structure). In this case the dynamics is determined by 
the bridge tension between the globules and the hydrodynamic friction experienced by the globules. 

We should stress that the whole consideration ignores up to now two important things: (i) hydrodynamical interac- 
tion and (ii) topological (or entanglements) effects. The hydrodynamical interaction can be simply taken into account 
by treating the solvent as a incompressible Navier - Stokes liquid coupled with chains monomers |2S|] . Within the 
self - consistent Hartree approximation the effective friction coefficient depends from the monomer - to - monomer 
correlation function, Q(n,t), so that the closed equation of motion for the Rouse modes correlator can be solved 
numerically. This point will be discussed in detail in a forthcoming publication. 

The systematic inclusion of topological constraints in the dynamics of collapse (or swelling) is a much more com- 
plicate theoretical problem. It was argued in ref. that at the chain length N ^> N c (where N c is an effective 
entanglement length) the topological constraints become essential. The self similar process of crumpling which we 
have discussed in Sec.IIA is restricted also by topological constraints as soon as a collapsed segment is longer then 
N c . These segments have a fixed topology and in a sense are in a partially equilibrium state called "fractal crumpled 
globule". The characteristic time of of this first stage of collapse is £ C oiiapsc ~ ^coiia P sc(l + b 6 /N c w) , where w is a 

third virial coefficient and Collapse ^ s *h e collapse time without the topological effects. At the subsequent stage the 
"crumpled globule" relaxes to the full equilibrium via the penetration of the chain ends through the globule and many 
knots formation. This stage has a reptational mechanism and as a result the characteristic relaxation time is scaled 
as itop ~ A 3 . Unfortunately, nowadays it is not quite clear how to incorporate the topological constraints in the 
equation of motion, so that this is a challenging problem for the future discussion. Finally we can mention as an 
interesting possible application of our approach the problem of the globule dissolution under an external force |?7j . 
This problem is also important for the interpretation of the stress - strain relations in polymer networks p8[ . 
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In order to accomplish the RPA - calculation let us make as usual |2J] the transformation to the collective solvent 
density 
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APPENDIX A: INTEGRATION OVER SOLVENT VARIABLES 
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(Al) 



p=i 



and response field density 
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d 




(A2) 



p=i j=i 
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These transform the influence functional (3.6) to the form 



R, R 



= In / Dp(k,t)Dn(k,t) 



t)V ss (k) 



{f f d d k 
W[p,n]- J dt J ^-p7r(-k,<)p(k, 

f f d d k 

+ Z2J dtiR j( s ^) J j^yikjV ps (k)p{k,t)ex V {ikR(s,t)} 



s=l ' 
N 



f f d d k 

£j dt J ^p7r(k,t)7 ps (-k)exp{ikR(s,t)} 



(A3) 



where the functional W depends only on properties of the free system 

M 



l(0) 
olvent 



W{p,n} = In / Y[ Dr ip) Dr{p) exp |a 
J P =i 

X 5^p(T,t)-J^6(T-TM(t)fj 

' M d 

7r M)-££ irf (t)V^(r-rW(t) 



} 



x 5 



p=i j=i 



(A4) 



Here A 



(0) 

solvent 



is the free solvent action. Following ref. pj 



one can expand W{p, 7r} up to the second order 



with respect p and tt, which formally corresponds to the dynamical RPA. Then the solvent variables in (A3) can be 
integrated over and for GF we obtain the following result 



Z{---} = J DR J {s,t)DR J (s,t) expjA cff R(s, t), R(s, i)j | 



(A5) 



where 



leff 



R, R 



- 4(°) 
1 

+ 2 



1 JV 



d d k 

W) 



2 k 3 k p \V vs {k)\ 2 e^{ik [R(l)-R(l')]}^(l')Soo(k,t-i') 



zAj(l) 



d d k 



ikjv(k) exp {ik [R(l) - R(l')]} 6(t - t') 



1 JV JV 



(2ir) d 3 

d d kd d q . 



(27T) 



ifcjii;(k, q) 



2d '""J 



exp {zk [R(l) - R(l')] + iq [R(l') - R(l")]} - *W - *") . 



(A6) 



where we have used a short hand notations: J 1 = Yl^—i dt, 1 = (s, t) and -Aj^ ra is the free chain action. In 
eq.( A6) )5oo(k, t) and Soi(k, t) stands for the solvent RPA - correlation and response functions correspondingly. They 
take especially simple forms after Fourier transformation with respect to time argument. This can be written as 



Sbo(k,u;) 



$ 00 (k,w) 



and 



S Q i(k,uj 



[1 + V M (k)S 10 (k,w)] [1 + F ss (k)$ i(k, to) 
$oi(k,w) 



[l + y ss (k)$ i(k,w) 



(A7) 



(A8) 
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where 3>oo(k, u>), $01 (k, u>) = ^^(k, ui) are correlation and response function for the free solvent correspondingly. 
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